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There is experimental evidence for a direct CP asymmetry in singly Cabibbo suppressed D decays, 
AAcp ~ 0.006. Naive expectations are that the Standard Model contribution to AAcp is an order 
of magnitude smaller. We explore the possibility that a major part of the asymmetry comes from 
supersymmetric contributions. The leading candidates are models where the flavor structure of the 
trilinear scalar couplings is related to the structure of the Yukawa couplings via approximate flavor 
symmetries, particularly (7(1), [(7(1)] 2 and U(2). The recent hints for a lightest neutral Higgs boson 
with mass around 125 GeV support the requisite order one trilinear terms. The typical value of 
the supersymmetric contribution to the asymmetry is AA^p^ ~ 0.001, but it could be accidentally 
enhanced by order one coefficients. 



INTRODUCTION 



THE SUPERSYMMETRIC PARAMETERS 
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The world average for the direct CP asymmetry in 
singly Cabibbo suppressed D decays, based on measure- 
ments by E687 [3L CLEO RB, E791 Q, FOCUS @, 
BaBar Q, Belle 0, CDF g g] and LHCb [13], is now 
4.3(7 away from zero [Til ] : 



= -0.00656 ±0.00154. 



Here. 
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In AAcp, that is the difference between asymmetries, 
effects of indirect CP violation largely cancel out [l2j ]. 
Thus, AAcp is a manifestation of CP violation in decay. 

The Standard Model (SM) contribution to the individ- 
ual asymmetries is suppressed by a CKM factor of order 

21m ( ^-yi ^ ~ 1-2 x 10~ 3 , and by a loop factor of order 

a s (m c )/Tr ~ 0.1. While one cannot exclude an enhance- 
ment factor of order 30 from hadronic physics [l3l - T2l| , in 
which case (p} will be accounted for by SM physics, it 
is interesting to explore the possibility that new physics 
contributes a major part of AAcp- 

The size of new physics contributions to AAcp is often 
constrained by other flavor-related observables, such as 
mixing or e'/e [22j. Supersymmetric models, via 



their contribution to the chromomagnetic operator, can 
generate large enough asymmetry in D decays without 



The 6x6 mass-squared matrix for the up- and down- 
type squarks can be decomposed into 3x3 blocks, q — 
u, d, 



\Ml\Ml\ 



D, F-terms, 



(3) 



where L and R denote SU(2) doublets and singlets, re- 
spectively. We denote the average squark mass by m. 
Then, it is convenient to parameterize the supersymmet- 
ric contributions to flavor changing processes in terms of 
dimensionless parameters, 



(4) 



where M, N = L,R. When, to a good approximation, 
only two squark generations are involved, one can express 
these parameters in terms of the supersymmetric mixing 
angles, (K^ 4 )ij, and the mass-squared splittings between 
squarks, Am 2 ^: 



Am 2 



(5) 



The parameters that are most relevant to AAqp are 
&ll = {$l L )i2 and S lr = {Sl R )i2, wni ch generate the 
chromomagnetic operator with Wilson coefficient given 
by 



C Sg = F(x)6ll + G(x)^5 lr , 



(6) 



conflicting with these observables 
work, we investigate whether this scenario is likely to 
be realized in supersymmetric models with viable and 
natural flavor structure. 



12|, |23j, |24|. In this where x = {mA/fh ), and the functions F and G can be 



found, for example, in Ref. (l2j. Given that G(x) is larger 
than F(x) by a factor of a few, and the enhancement fac- 
tor of wig/nic, the dominant contribution in the models 
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that we consider comes from 5lr- It can be estimated as 
follows 23[ : 



AA: 



SUSY 
CP 



0.006 



lm(S LR ) ITeV 



(7) 



0.001 m 
In the following sections, we investigate whether 

Im{5 LR ) ~ 0.001 (8) 
can plausibly arise in supersymmetric flavor models. 

SUPERSYMMETRIC FLAVOR MODELS 

If the soft supersymmetry breaking terms had a generic 
flavor structure ("anarchy"), the supersymmetric con- 
tributions to flavor changing neutral current processes 
would exceed experimental constraints by orders of mag- 
nitude. Thus, these terms must have a special structure. 

The most extreme solution to this "supersymmetric 
flavor puzzle" is a constrained version of minimal flavor 
violation (MFV): at the supersymmetry breaking medi- 
ation scale, squark masses are universal and the trilinear 
scalar couplings (the A terms) are proportional to the 
corresponding Yukawa matrices. Such a situation arises 
naturally in various mediation schemes, most notably 
gauge- and anomaly-mediated supersymmetry breaking. 
The renormalization group evolution does generate fla- 
vor changing effects in the soft breaking terms, but MFV 
implies a very strong flavor suppression of these effects, 



$lr oc — (V us V* s y s 
m 



V ub V; b y 2 b ) < O(l0- 7 ). (9) 



(For the exact expression in anomaly mediation, see [25j.) 

It is possible, however, that the supersymmetric flavor 
structure is related to that of the Standard Model, but is 
not MFV. This is the case in models where an approxi- 
mate flavor symmetry dictates the structure of all flavor 
changing couplings. In what follows, we examine several 
such symmetries - 17(1), [[/(l)] 2 , U{2) and [U{2)f - with 
regard to their implications for Xm(SLR)- 



Abelian Symmetries 

The Froggatt-Nielsen (FN) framework (2|| postulates 
an approximate U(l) symmetry, broken by a spurion with 
value <C 1. Assigning different charges to the different 
quark generations results in parameterically suppressed 
quark mass ratios_ and mixing angles. In supersymmet- 
ric FN models 
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28j , the squark spectrum is anarchi- 
cal, up to some level of degeneracy between the first two 
generations from renormalization group evolution (RGE) 
effects, but the mixing angles are small. 

With a single U(l) symmetry, the parametric suppres- 
sion of squark flavor parameters is related to quark flavor 



parameters, independent of details of the model such as 
the size of the spurion and the charge assignments. In 
particular, the following relations hold for the entries that 
are relevant to c — > u transitions [29j : 



( S LL)l2 
(<^?fl)l2 



\Vus 
T3 



rzm c \V us y 

a m c \Vus\ 
in rh 



(10) 
(11) 
(12) 



where a is the typical scale of the A-terms. The l/r% 
factor, defined in Ref. [13], represents the gluino- related 
RGE effect which generates some level of degeneracy be- 
tween the first two squark generations, Am 2 2 /m 2 ~ I/V3. 

Throughout this paper we assume unsuppressed CP 
phases and x ~ 1. The supersymmetric contribution to 
ck that is proportional to (<5^ L ) 12 (S RR ) 12 ~ rrid/ (r§7n s ) 
is too large, unless > 440 (TeV/m) (or, equivalently, in 
the language used by CMSSM practitioners, m 1 / 2 /mo > 
7). Assuming that this is indeed the case, the model is 
viable, and provides 



Im{5 LR ) ~ 1.5 x 10" 



a 1 TeV 



(13) 



Comparing Eq. (fT3"|) to Eq. © , we learn that for super- 
symmetry to account for AAqp, the ratio (a/rh) should 
be large. Taking ao, rriQ and TO1/2 to stand for the Planck 
scale values of the j4-terms, squark masses and gluino 
mass, respectively, and using the approximations of Ref. 
[301, we obtain 



3a 



rh m ^/l + 8{m 1/2 /m ) 2 



ao/mx/2 (mi/2 > m ), 
3a /m (m 1/2 < m ). 

(14) 



> 



Given that the Z7 (1) models are only viable if mi/ 2 
7mo, the optimal enhancement occurs for ao > mi/ 2 3> 
to which might, however, lead to negative squark 
masses-squared. 

The single U(l) models lead to the following simple 
parametric relation between the up and down sectors: 



(S d LR )i2 m s 



(15) 



The (<5f/j)i2 parameter is constrained, however, by e'/e: 
(Kr)i2 ^ 4 x 1°" 5 (m/TeV) (see [HIHl and references 
therein). We thus obtain an upper bound that is inde- 
pendent of the flavor-diagonal scales, 



(Slr)™ <5x 10- 



TcV 



(16) 



Moreover, the approximate symmetry relates (S^ R )i2 
to flavor diagonal parameters, 



^ TO (^Efl)i2 ^ m c \V us \ 

Z m ( S LR)ll m q 



(q = u,d). 



(17) 
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Assuming phases of order one (which we must do to 
explain AAcp), these flavor diagonal parameters are 
bounded by electric dipole moment (EDM) constraints 



(see [3(1 and references therein), (S^ R )u < 3 x 

- 6 (m/TeV). The 



(18) 



10~ 6 (m/TeV) and (Sf R )u < 2 x 10" 
resulting bounds are 

(^)i2<3xl0" 4 ^ (from (SIM, 

(^)i2<8xl0- 5 ^ (from (On)- (19) 

We conclude that FN models with a single U(l) are un- 
likely to account for AAcp 3> 0.001. Of course, since the 
FN mechanism only dictates the parametric suppression, 
it is impossible to exclude an accidental enhancement of 
(^E_r)i2 by the order-one coefficient. 

Models with an FN symmetry ][/(l)] 2 allow one to take 
advantage of the holomorphicity of the superpotential 
to obtain vanishing entries in the Yukawa and A matri- 
ces and to strongly suppress entries in the squark mass- 
squared matrices (compared to the single U (1) case). 
This feature was first employed in Refs. [2jJ, |28| to align 
in a very precise way the squark and quark mass matrices 
in the down sector. 

The flavor structure of the Yukawa and A terms in 
these models can be written as follows: 



Y° 



av d 



M 2 LR 
av u 



Vd Vb\V ub \ ^ 

y s Vb\V cb \ , (20) 
y b j 

Vu Vc\V U a\ Vt\V u b\ 

*2 Vc Vt \V cb \ I . (21) 
Y&y c /\V A \ y t 



The four holomorphic zeros in the down sector are es- 
sential to obtain an effective alignment 29]. The (21) 
and (31) entries in the up sector either (i) get their naive 
parametric suppression (of order y u /\V U8 \ and y u /\V u b\, 
respectively) or (ii) vanish. 

In both cases, the contribution to €k from 
(<5 'll) 12 ($rr) 12 for unsuppressed phases is too large, un- 
less RGE generates degeneracy, r$ > 18 (TeV/m). (For 

moderate suppression of phases, D° — Z?° mixing pro- 
vides the strongest constraint, requiring a milder degen- 
eracy 32]: In case (i), the estimates (fTU)) and ([TT]) hold, 
and the contribution from (S^ L )i2(S RR )i2 is too large, 
unless RGE generates degeneracy, r 3 > 7 (TeV/ra) [301 ] . 
In case (ii), {S RR )i2 is suppressed compared to Eq. (fTTj) . 
and the contribution from ](<5ii)i2] 2 only requires a very 
mild degeneracy.) 

In either case, the parametric suppression of (S^ R )i2 
is as in Eq. (|12[) , and the numerical estimate is as in Eq. 
(|13p . The parametric relation between the up and down 
sectors of Eq. (fT5|) does not hold in the [£/(!)] 2 models 



since (<^^)i2 is further suppressed, and so the constraint 
of Eq. ((T6| does not hold. The constraints of Eqs. ([18]) 
and (|19|) hold, leading to the same conclusion as in the 
single U(l) case: The parametric suppression is such that 
the contribution to AAqp falls an order of magnitude 
short compared to the benchmark value of Eq. ©. How- 
ever, one cannot exclude the possibility that: (i) the large 
hadronic uncertainties in the EDM calculation are such 
that the bounds are weaker by an order of magnitude; 
or (ii) the order one uncertainty from the unknown coef- 
ficients provides further enhancement; or a combination 
of the above. 



Non-Abelian Symmetries 

In U (2) models, the first two generations are in a dou- 
blet and the third generation in a singlet of the U(2) 
symmetry [33j-|36|. With a two-stage symmetry break- 
ing, the structure of the Yukawa, A and M 2 matrices is 
as follows (see, for example, [361]): 



M 



2q 

LB 



M 



2q 
NN 




m\(l + el) e 2 mf 



(22) 



V o 



As in the Abelian case, all non-vanishing entries have 
unknown coefficients of order one, but Y^ 2 = —Y 2 \ and 
there are relations for the M^ N matrices that follow from 
hermiticity. 

As concerns the S^ R parameters, their parametric sup- 
pression is similar to the £7(1) model. Hence, Eqs. (|12l) . 
(HH), (HHJ) and (HU) all hold. 

The main phenomenological difference of the U(2) 
model with respect to the U(l) model is that the first 
two squark generations are quasi-degenerate already at 
the mediation scale, with a mass splitting Ato 2 2 /to 2 ~ 
e 2 ~ 10~ 3 . Hence, the model is viable even without in- 
voking flavor-universal RGE effects. 

A flavor ]C/(2)] 3 symmetry [37f is motivated by the ten- 
sion between the measured value of the CP asymmetry 
S^k, and its theoretical value in the Standard Model 
derived from a global CKM fit. To alleviate this ten- 
sion, a new physics contribution to B° — B mixing of 
order 10 percent of the total amplitude is required. In 
a U (2) model, such a contribution entails a contribution 
to tK of order 100 percent, which is unacceptable. A 
U(2)q x U(2)u x U(2)d model, with minimal spurion 
content - 7(2,1,1), AY„(2,2,1) and AY d (2,l,2) - al- 
lows one to suppress the contribution to €k ■ 
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The structure of the Yukawa and .4 matrices is as fol- 
lows (37|: 

where y q3 — yt(yb) for q = u(d), x q is a complex free 
parameter of order one, AY q is a 2 x 2 matrix, and V is 
a 2 x 1 vector. 

This structure is quite unique in that one and the same 
spurion, AY q , determines the structure of the 2x2 up- 
per left block of both Y q and M L q R . Consequently, to 
leading order in the breaking parameters, {5 q LR )i2 = 0, 
and the supersymmetric contribution to AAqp vanishes. 
Corrections to Slr arise at the order y c ytV* b V u b and are 
negligible. 



close to the upper bound from EDMs, Eq. (p~9|) . We 
note that it is possible in specific U(l) 2 models to further 
suppress the contribution to the kaon system, such that 

the strongest bound comes from D° — D mixing and 
gives r/r 3 < 0.8 (m/TeV), relaxing the constraint (|27|) 
by a factor of 4. 

Both £/( 2 ) and [U{2)] 3 models do not require further 
flavor suppression, and are viable for r < 1, with predic- 
tions as in the non-hybrid models. 

We conclude that hybrid models with a [C/(l)] 2 , U(2) 
or U(2) 3 symmetry generate AAqp of the same size as 
non-hybrid models. The size of 5lr allowed by hybrid 
models with a single U(l) is somewhat smaller. 

A TERMS AND THE LIGHTEST HIGGS MASS 



Hybrid Mediation 

In models of hybrid mediation, the dominant source 
of supersymmetry breaking is MFV, but there are non- 
negligible contributions from Planck scale physics that 
do not obey the MFV principle. Examples include high- 
scale gauge mediation [38| and a class of models with 
anomaly mediation [39[ . At the messenger scale, the rel- 
ative size between the soft masses-squared arising from 
gravity and MFV physics is given by r ~ TO g r av/ m MFV- 



The c — > u couplings are given by 



(6 u LL )i2 

i S RR)l2 



r\V U s\ 

1 

r to, 



r3m c \V us \ ' 



(24) 
(25) 



and the expressions for 5lr remain as in the pure gravity 
case, Eqs. t$TZ§ and (fT3|) . 

The EDM constraints of Eqs. ([18]) and ([19]) hold. One 
can now ask what further constraints arise when linking 
the trilinear terms and the soft masses as is characteristic 
in hybrid models. If all gravity soft terms are dictated 
by a single scale, then clq ~ y/r tomfVj where tomfv is 
the typical messenger scale MFV soft mass. The relevant 
combination entering 6lr is then a/rh < 3yr/r3, where 
the numerical factor stems from RGE and is largest for 
high scale mediation. In the following we analyze this 
single gravity scale scenario in the FN context. 

In models with a single U(l) and order one phases kaon 
mixing requires r/r^ < 0.002 (rh/TeV) [3(j. Therefore, 



lm((5 M ) < 0.2 x 10" 



TeV 



(26) 



a stronger constraint than the EDM bound of Eq. |Tt 

In [C/(l)] 2 models the kaon system constrains r/r 3 < 
0.06 (m/TeV), and so 



Zm(6 LR ) < 1 x 10" 



TeV 



(27) 



In supersymmetry AAqp can be interpreted via the 
left-right mixing Slr which requires unsuppressed trilin- 
ear couplings with respect to the squark masses, a/rh ~ 
0(1), see Eq. (jT5]) . At the same time the recent hints 
from ATLAS and CMS of a neutral Higgs boson with 
mass near 125 GeV [40( implies that the stops - if not 
decoupled - are largely mixed as well [4l|, l42j : 



\At/vt-n/taa0\>M s , 



(28) 



where Ms denotes the geometric mean of the stop 
masses. In the FN models, where the flavor structure 
of the A terms is parametrically similar to that of the 
Yukawas, A ~ Y, the stop A terms at the weak scale can 
be written as [4311 



At/Vt 



Aa, Aa = y\ba^ 



H/2C, 



(29) 



with positive RGE-induced coefficients b, c of order one. 

For positive fi or sufficient tan /3 suppression one ob- 
tains from Eq. (|28j) a lower bound that supports a size- 
able supersymmetric AAcp, 



a I m > Ms I to 



(30) 



for A t > and unsplit spectrum where the stops are not 
too far away from the other squarks. Negative A t < can 
arise in scenarios with tiny or vanishing eto such as gauge 
mediation, which lead to acceptable phenomenology only 
for sufficiently large gluino masses. 

While the Higgs signal needs to be consolidated, it 
is interesting that if confirmed, the current mass of 
~ 125 GeV points to a similar region in supersymmetric 
parameter space as the interpretation of AAcp. 



CONCLUSIONS 

Supersymmetric models can contribute to direct CP 
violation in singly Cabibbo suppressed D decays at the 



5 



level observed by experiments, AA^, SY ~ 0.006, without 
conflicting with phenomenological constraints from D° — 
D mixing or e'/e. This is naturally the case if the flavor 
changing parameter (S LR )i2, generated by trilinear scalar 
couplings, is of order 10~ 3 . 

In minimally flavor violating supersymmetric models, 
such as those of gauge mediation and anomaly mediation, 
(<^l_r)i2 is orders of magnitude too small. Thus, to ac- 
count for A Aq p , one has to go beyond minimal flavor vi- 
olation. We examined models where the flavor structure 
of the soft breaking terms is dictated by an approximate 
flavor symmetry. 

We found that quite generically in such models, 
(^Lfi)i2 i s flavor-suppressed by (m c | V us \/rh), which is of 
order a few times 10 -4 . There is however additional de- 
pendence on the ratio between flavor-diagonal parame- 
ters, a/ in, and on unknown coefficients of order one, that 
can provide enhancement by a factor of a few. 

In most such models, however, the selection rules that 
set the flavor structure of the soft breaking terms, relate 
0^)12 to (5 d LR )i2 and to (5\ 
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5u,d\ 



which are bounded 



from above by, respectively, e'/e and EDM constraints. 
Since both e'/e and EDMs suffer from hadronic uncer- 
tainties, small enhancement due to the flavor-diagonal 
supersymmetric parameters cannot be ruled out. Addi- 
tionally, it is still possible that (6^ R )i2 is accidentally 
enhanced by the order one coefficient. 

Chirality-flipping couplings between the first and sec- 
ond generation up squarks can effectively arise also via 
chirality-flipping in the third generation 23[. In all fla- 
vor models considered here, the effective Slr generated 
is at most parametrically of the same order as the direct 
contribution, accompanied by an additional l/r|. For 
t/(l) and [i7(l)] 2 models this provides extra suppression, 
while for U(2) models the effective contribution to Slr is 
flavor-suppressed with respect to the direct one. In any 
event, the constraints remain the same and the analysis 
stands. 

We conclude that it is possible to accommodate 
AAqp ~ 0.006 in supersymmetric models that are non- 
minimally flavor violating, but - barring hadronic en- 
hancements in charm decays - it takes a fortuitous ac- 
cident to lift the supersymmetric contribution above the 
permil level. 
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